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ABSTRACT 

The steady-state structure of self-gravitating, magnetized accretion disks is studied using a set of 
self-similar solutions which are appropriate in the outer regions. The disk is assumed to be isothermal 
and the magnetic field outside of the disk is treated in a phenomenological way. However, the internal 
field is determined self-consistently. The behaviour of the solutions are investigated by changing the 
input parameters of the model, i.e. mass accretion rate, coefficients of viscosity and resistivity, and the 
magnetic field configuration. 

Subject headings: accretion, accretion discs - hydrodynamics 



1. INTRODUCTION 

One of the key physical ingredients in accretion disks is 
self-gravity. It plays a significant role in many such sys- 
tems, ranging from protostellar disks to Active Galactic 
Nuclei (AGN). The radial and vertical equations for the 
disk structure are significantly modified because of the pos- 
sible impact of the self-gravity, although traditional mod- 
els of accretion disks ignore self-gravity just for simplicity 
(e.g., Pringle 1981). Nevertheless, Deviations from Kep- 
lerian rotation in some AGN and the flat infrared spec- 
trum of some T Tau stars can both be described by self- 
gravitating disk models. 

Because of the complexity of the equations and in order 
to obtain analytical results, authors have studied the ef- 
fects related to disk self-gravity either in the vertical struc- 
ture of the disk (e.g., Bardou et al. 1998) or in the radial 
direction (e.g., Bodo & Curir 1992). In a situation where 
we lack strong empirical evidence for the detailed mecha- 
nisms involved, Bertin (1997) proposed a new class of self- 
gravitating disks for which efficient cooling mechanisms 
are assumed to operate so that the disk is self-regulated 
at a condition of approximate marginal Jeans stability. 
Thus, the energy equation is replaced by a self-regulation 
prescription. He showed that in the absence of a central 
point mass, there is a set of self-similar solutions describing 
the steady-state structure of such self-regulated accretion 
disks. The self-similar solution corresponds to a flat ro- 
tation curve, while the disk has a fixed opening angle. In 
fact, Bertin's self-similar solution is a generalization of the 
self-similar solution for self-gravitating disks (Mestel 1963) 
dominated by viscosity. 

Subsequent analysis confirmed the validity of this sim- 
ple self-similar solution (Bertin & Lodato 1999). More- 
over, such kinds of self-regulated accretion disk can suc- 
cessfully describe the spectral energy distribution of proto- 
stellar disks (Lodato & Bertin 2001). Recently, Lodato & 
Rice (2004) studied the transport associated with gravita- 
tional instabilities in a relatively cold disk using numerical 
simulations. They showed that the disk truly settles into a 
self-regulated state, where the axisymmetric stability pa- 
rameter Q « 1 and where transport and energy dissipa- 
tion are dominated by self-gravity. However, these studies 
neglected the possible effects of magnetic fields. Many au- 



thors tried to construct models for magnetized disk either 
in a phenomenological way based on some physical con- 
siderations (e.g., Schwartzman 1971; Lovelace et al. 1994; 
Kaburaki 2000; Shadmehri 2004; Shadmehri & Khajen- 
abi 2005) or by direct numerical simulations (Stone et al. 
1999; McKinney & Gammie 2002). Bisnovatyi-Kogan & 
Lovelace (2000) suggested that recent papers discussing 
Advection-Dominated Accretion Flow (ADAF) as a possi- 
ble solution for astrophysical accretion should be treated 
with caution, particularly because of our ignorance sur- 
rounding the magnetic field. Models of magnetized ac- 
cretion disks with externally imposed large scale vertical 
magnetic field and anomalous magnetic field due to en- 
hanced turbulent diffusion have also been studied (see, 
e.g., Campbell 2000; Ogilvie & Livio 2001). These models 
are restricted to subsonic turbulence in the disk and the 
viscosity and magnetic diffusivity are due to hydrodynamic 
turbulence. 

Recently, we presented a set of analytical self-similar 
solutions for the steady-state structure of a magnetized, 
radiation-dominated disk (Shadmehri & Khajenabi 2005). 
Although this analysis is just a simple model, one can see 
some possible effects of the magnetic field on the structure 
of radiation-dominated disks, at least at a fundamental 
level. However, we applied some simplifying assumptions 
concerning the magnetic field both outside and inside of 
the disk. Our approach was based on interesting studies 
by Lovelace et al. (1987) and Lovelace et al. (1994). We 
noticed that it is possible to take a similar approach to 
study magnetized, self-gravitating disk. 

In this paper, we relax the self-regulation condition of 
Bertin (1997) simply by assuming an isothermal equation 
of state. We present self-similar solution for the steady- 
state structure of such self-gravitating disk for which 
the effect of the magnetic field is taken into account. 
While the magnetic field outside of the disk is studied 
in a phenomenological approach similar to Lovelace et al. 
(1994), the magnetic field inside the disk is obtained self- 
consistently. We show that magnetic fields may cause sig- 
nificant changes in the typical behaviours of the solutions 
compared to the nonmagnetized case. The basic equations 
of the model are described in section 2. Self-similar solu- 
tions are obtained in section 3. In spite of the similarity 
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of Bertin's solutions (1997) to ours, there are differences 
between his work and the present analysis. A summary of 
results are discussed in the final section. 

2. GENERAL FORMULATION 

The basic equations arc integrated over the vertical 
thickness of the disk. The mass continuity equation be- 
comes 

-27ri?SwR = M, (1) 

where vn is the radial velocity and E = J dzp ~ 2/ip is 
the surface density of the disk. The half-thickness is de- 
noted by h, where we consider the magnetic field effect on 
the disk thickness. We will see that the disk can be com- 
pressed or flattened depending on the field configuration. 
Also, wc note that since the radial velocity is negative for 
accretion (i.e., wr < 0), the accretion rate M as an input 
parameter of our model is positive. 
The radial momentum equation reads 



^ dv^ vl _ dP ^ GMjR) 



Fr^dz, (2) 



where P = J dzp is the integrated disk pressure and M{R) 
is the mass of the disk within radius R . Wc are neglect- 
ing the mass of the central object comparing to the disk 
mass. This is relevant to protostellar disks at the begin- 
ning of the accretion phase during which the mass of the 
central object is small and self-gravity of the disk plays 
an important role. Also, one may argue that our model 
corresponds to disks at large radii because the effects of 
the central mass become unimportant in the outer regions 
of the disk. The term associated with self-gravity of the 
disk, i.e. GM{R) / R^ , is generally applicable to spherical 
geometry. In fact, the correct term takes on an integral 
form (see equation (4) of Bertin 1997). But, as we will 
show, our solution describes a disk in which the surface 
density is inversely proportional to radius. Interestingly, 
for such a disk which resembles a Mestel disk, the cor- 
rect gravitational force is known which is exactly equal 
to our approximate formula (Binney & Tremaine 1987). 
It implies that we consider the exact gravity force in our 
equations. Now, we can write 



dM{R) 
dR 



(3) 



The last term of equation (2) represents the height- 
integrated radial magnetic force which can be written as 
(Lovelace et al. 1994; hereafter LRN) 

„ . Id 

^ [h<Bl>] + i-^ {Bl + Bl - Bl\, (4) 
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where < •■• >= J_j^dz{- ■ ■)/{2h), and the h subscript 
denotes that the quantity is evaluated at the upper disk 
plane, i.e. z = h. Similarly, integration over z of the 
azimuthal equation of motion gives 

(5) 



where 

R^F:;^^Hz = ^{R^B^B,)^ - l.^{R^BnB^)^ 

Here, the last term of equation (5) represents the height- 
integrated toroidal component of magnetic force multiplied 
byi?2. 

In our model, we also assume (Shakura & Sunyaev 1973) 

I' = aCsh, (7) 

where Cg is the local sound speed and a is a constant less 
than unity. 

The z component of equation of motion gives the condi- 
tion for vertical hydrostatic balance, which can be written 

as 

= ttGS^ + i^)[{Bn)l + [BX] - ^ (^H)h^. (8) 

where the term on the left hand side of the above equa- 
tion corresponds to the thermal pressure. Since we neglect 
the central object, the self-gravity of the disk in the ver- 
tical direction leads to the first term on the right hand 
side of the above equation. Now, we can treat the internal 
magnetic field using the induction equation. LRN showed 
that the variation of B^, with z within the disk is negligible 
for even field symmetry. Moreover, B^ and B^ are odd 
functions of z and consequently dBp^/dz « {BYi)h/h and 
dB^/dz « (B^)hlh. Krasnopolsky & Konigl (2002) ap- 
plied similar configurations to study the time-dependent 
collapse of magnetized, self-gravitating disks. Thus, 

B^{R,z) = ^{Bn)h,B^{R,z) = ^{B^, 

B,{R,z)=B,{R), (9) 

and the induction equation reads 

-i?i?,.;R - ^(SR)h + = 0, (10) 

where the magnetic diffusivity 77 has the same units as 
kinematic viscosity. We assume that the magnitude of r] 
is comparable to that of the turbulent viscosity u (e.g., 
Bisnovatyi-Kogan & Ruzmaikin 1976; Shadmehri 2004). 
Exactly in analogy to the alpha prescription for v, we are 
using a similar form for the magnetic diffusivity 77, 

T] = T]oCsh, (11) 

where ijo is a constant. Note that rj is not constant and 
depends on the physical variables of the fiow, and in our 
self-similar solutions, as we will show, rj scales with radius 
as a power law. This form of scaling for diffusivity has 
been widely used by many authors (e.g., Lovelace, Wang 
& Sulkanen 1987; Lovelace, Romanova & Newman 1994; 
Ogilvie & Livio 2001; Riidiger & Shalybkov 2002). 

While equation (10) describes the transport of a large- 
scale magnetic field (here, Bz{R)), the values of (i?R)h and 
{B^)h are determined by the field solutions external to the 
disk. Instead, we are following the approach of LNR, in 
which the external field solutions obey the relations 

(-BrJii = f^iBz, (-By)h = l^ipBz, (12) 
where /3r and are constants of order unity {(3^ < 0). 
Thus, one can simply show that < B^ >= (3^B'^/3, 
< Bl >= PlB^/3 and < B^B^ >= (irP^B^/3. Note 
that these beta factors are really empirical scalings which 
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now have abundant support from computer simulations 
of MHD and Poynting outflows from disks. Ustyugova et 
al. (1999) gives a detailed analysis and provides strong 
evidence for /J,. ^ 1/3^1 ~ 1- 

To close the equations of our model, we can write 
the self-regulation condition as has been done by Bertin 
(1997). But one should naturally expect another form of 
self-regulation prescription in the magnetized case. In fact, 
self-regulation results from a competition of dissipation 
and instabilities. Indeed, the instability of a magnetized 
disk is different from that of unmagnetized disk. Effects 
of the magnetic field on linear gravitational instabilities 
in two-dimensional differentially rotating disks have been 
investigated in detail by Elmegrecn (1987, 1994), Gam- 
mie (1996) and Fan & Lou (1997). Axisymmetric and 
nonaxisymmetric perturbations show significantly differ- 
ent behaviours. While axisymmetric instability in a thin 
disk borrows its physical groimds from Toomrc (1964) 
and requires Q < 1, the literature has lacked an explicit 
theoretical evaluation of critical Q for nonaxisymmetric 
gravitational runaways. In order to avoid such difficul- 
ties one can consider the axisymmetric stability. The 
presence of the magnetic fields modifies the Toomre cri- 
terion for a hydrodynamic disk in such a way that mag- 
netized disks are unstable to axisymmetric perturbation if 
Qm = Q(l + 1//3) < 1, where /3 is the ratio of the thermal 
pressure to the magnetic pressure (e.g., Shu 1992). 

There is a fundanicaital point concerning the self- 
regulation prescription. When the disk mass becomes 
large enough to induce global instability there is nothing 
like self-regulation taking place. We can see numerical sim- 
ulations which describe massive unstable disk (e.g., Bon- 
nell 1994; Matsumoto & Hanawa 2003). However, consid- 
ering existing numerical simulations, we think, it is not 
possible simply to say that all massive disks fragment be- 
cause, to our knowledge, not only current numerical simu- 
lations suffer from their own limitations, but also they do 
not generally give a fully consistent picture as for global 
fragmentation. This situation becomes more complicated 
when we consider magnetic fields. It seems that there is a 
complicated interaction between gravitational instability 
and MHD turbulence that infiuences disk structure, but 
MHD turbulence reduces the strength of the gravitational 
instability (Fromang 2005). 

Almost all numerical simulations show that global in- 
stabilities are highly depend on the thermal state of the 
disk (Gammie 2001, Johnson &: Gammie 2003). However, 
thermal state of the disk has been neglected by some au- 
thors and so their approach is not a complete analysis. 
Simulations by Mastsumoto & Hanawa (2003) which do 
not include detailed thermal evolution, predict fragmen- 
tation in an early phase. But all these fragments are on 
tight orbits and are likely to merge due to disk accretion. 
Matzncr & Levin (2005) argue analytically that viscous 
heating and stellar irradiation quenches fragmentation and 
conclude that numerical simulations lead to fragmentation 
which do not account for irradiation and are unrealistic. 
On the other hand, not necessarily all numerical simula- 
tions of massive disks lead to fragmentation. For example, 
Pickett at al. (2000) failed to obtain fragmentation in 
similar thermodynamics condition. Thus, whether or not 
massive disks fragment due to global instability remains 



controversial. At this stage what we can say is that even 
if a massive disk with no star in the center does not frag- 
ment, it does not mean that self-regulation is occurring in 
the disk. In order to avoid such difficulties, we relax the 
self-regulation condition and replace it by an isothermal 
assumption, i.e. 



p = ^ci 



(13) 



where Cs is constant sound speed. We hope our simple ap- 
proach will be able to illustrate some possible effects of the 
magnetic fields on the structure of self-gravitating disks. 

3. SELF-SIMILAR SOLUTIONS 

Equations (1), (2), (3), (5), (8), (10) and (13) constitute 
the basic equations of our model for the steady-state struc- 
ture of self-regulated magnetized disk. We can solve these 
equations numerically using appropriate boundary condi- 
tions. However, before doing such an analysis it would be 
illustrative to study the typical behaviour of the solutions 
by applying semi-analytical methods, e.g. self-similarity. 
Indeed, any self-similar solution contains part of the be- 
havior of the system, in particular far from the boundaries. 
However the main goal of our analysis is just to illustrate 
the possible effects of the magnetic field on the steady- 
state structure of a self-gravitating disk. For this propose, 
we think, self-similar solutions are very useful. 

It is then straightforward to find a self-similar solution 
with the radial dependence 



Vo 

VRiR) 



= b, 



Vo 



PiR) 



^0 tlQ 



B.{R) 
Bo 

h{R) 

i?0 



Ho 



M{R) _ R_ 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



with a, b, c, rf, e, / and q being numerical constants which 
can be obtained from the equations (see below). Also, Rq, 
Eq, Pq, Vo and Bq are convenient units which reduce the 
equations into non-dimensional forms. We are assuming 

Vo = Cs = sJPo/^0 and Mo = Rl^o- 

If we substitute the above self-similar solutions in the 
main equations of the model, the following system of di- 
mensionless algebraic equations are obtained which are to 
be solved for o, 6, c, d, e, and /: 



oc = TO, 



-{ab^ + d) = -2a2 + [2(3, + -/(3 + fil - P^,)]e 



c)=/3^(2--//3,)e^ 



(21) 
(22) 

(23) 
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y =a2 + (//3r + /3r' + /3^)e^ (24) 

c-r?o/3r\/--?7o\/-/ = 0, (25) 
V a V a 

a, (26) 

where m, a, 7?0) A and are the input parameters. Also, 
m = M/{27tRo'EoVo) is the nondimensional mass accre- 
tion rate. The radial dependence of physical quantities 
both in the magnetized and the unmagnetizcd cases are 
similar. After some algebraic manipulations we can find 
a complicated algebraic equation for / and clearly only 
the real root is a physical solution. Having the real root of 
this equation, the other physical quantities are determined 
using the above equations. 

Now, we explore the parameter space of the input pa- 
rameters and their effects on the solutions. We restrict our 
study to values around unity for /S^ and 13^ (Ustyugova ct 
al. 1999). Since we find that the sensitivity of solutions on 
the parameter is not very strong, we plot all the phys- 
ical variables as a function of /3r for a fixed 0^. However, 
we comment about possible effects of variations of by 
the final figure. 

Figure 1 shows various physical variables as function of 
/3r with a = 0.1, rjo = 0.1, = —0.8 and the mass accre- 
tion rate is indicated on the plots. This figure helps us to 
see the dependence of the physical variables on the varia- 
tions of the mass accretion rate m. The plots on the top 
show surface density (left) and the opening angle (right) 
for m = 0.5, 1.0 and 2.0. We see as the parameter of the 
radial component of the magnetic field at the surface of 
the disk increases, the surface density decreases for a fixed 
mass accretion rate. But the surface density increases by 
increasing rh. Also, as the mass accretion rate increases, 
the disk becomes thinner. However, the disk thickness in- 
creases with /?r for a fixed accretion rate. The plots on the 
middle show typical behaviors of the radial and the rota- 
tional velocities. While the ratio of the radial velocity to 
the sound speed increases with increasing /3r, the ratio of 
the rotational velocity to the sound speed decreases keep- 
ing m constant. However, increasing the mass accretion 
rate causes the radial velocity to decrease. But for the 
rotational velocity, we see completely different behaviour, 
i.e. the ratio v,^/cs significantly increases with the accre- 
tion rate m. 

Our self-similar solution generally corresponds to > 1, 
where (3 is the ratio of thermal pressure to magnetic pres- 
sure at the surface of the disk. The plot on the bottom of 
Figure 1 (left) shows this behavior, although decreases 
with /Sj. for a fixed m. However, if the accretion rate in- 
creases the ratio of the thermal pressure to the magnetic 
pressure (3 increases which implies the effect of the mag- 
netic field on the structure becomes weaker. We have 
already discussed about the Toomre parameter Q and a 
modified version of this parameter, i.e. Qm because of 
the magnetic field. Figure 1 shows variation of both Q 
and Qm as function of /3r (curves corresponding to Qm are 
shown by dotted lines). Generally, we have Qm > Q be- 
cause of the dynamical effects of the magnetic field and as 
/3r increases these parameters increases to values closer to 
unity or even higher values. It implies a more stable disk 
according to the Toomre criteria as (3r increases. Also, the 



Toomre parameter rapidly increases when the mass accre- 
tion rate decreases. 

The r/o-dependence of the solutions is shown in the Fig- 
ure 2. In fact, this Figure is similar to Figure 1 except for 
changing the input parameter rjQ to different values 0.05, 
0.1 and 0.2. The other input parameter are the same as 
in Figure 1. The surface density S, the rotational veloc- 
ity and the ratio /3 decrease with increasing tjq. But 
the radial velocity vr, increases with increasing rjQ. We see 
that as 770 tends to the larger values, the disk thickness 
increases. Also, for a fixed tjo the disk thickness increases 
as Pr increases, however, the ratio h/R is more sensitive 
to the variation of P^- for higher values of rjQ. Interest- 
ingly, as the magnetic diffusivity coefficient rjo increases, 
the Toomre parameters Q and Qm significantly increase 
to values higher than unity. 

The a-dependence of the solutions is not strong as long 
as Pr w \p,^\ « 1. But we found that the ratio /? and the 
Toomre parameter change because of viscosity variations. 
Figure 3 (top, left) shows as a function of Pr- Clearly, as 
the viscosity parameter a increases, the ratio P increases. 
Also, the Toomre parameter increases with the viscosity 
parameter. We found that /3^-dependence of the solutions 
is weak for |/3^| « 1. However, there is some changes on 
the ratio p as P^, varies. Figure 3 (bottom, left) shows 
variation of /3 as a function of Pr for = —0.8, —1.0 and 
—1.2. Evidently, this ratio increases with |/3^|. Also, by 
increasing \p^\ the Toomre parameter increases, though it 
is not very significant. 

Bertin (1997) showed that the opening angle of self- 
regulated nonmagnetized disk depends only on Q, inde- 
pendent of the viscosity coefficient and the mass accretion 
rate. He argued that significantly large values of Q would 
be undesirable because they conflict with the thin disk 
approximation. But the disk thickness in our magnetized 
case depends on the mass accretion rate and the other in- 
put parameters as long as the beta factors are of order 
unity. For larger values of Pr or P^, the disk thickness 
may change because of the magnetic fields. 

4. CONCLUSIONS 

We have obtained a self-similar solution for a self- 
gravitating, magnetized, viscous disk. The solution has 
constant rotational and accretion velocities, independent 
of the radial distance. Since we are considering an isother- 
mal magnetized disk, the sound speed is constant. 

Our goal is to study the possible effects of the magnetic 
fields on the steady state structure of self-gravitating, mag- 
netized disk, at least at the physical level. Although we 
have made some simplifying assumptions in order to treat 
the problem analytically, our self-similar solution shows 
magnetic fields can really change typical behaviour of the 
physical quantities of a self-gravitating disk. Not only the 
surface density of the disk changes, but also the rotational 
and the radial velocities significantly change because of 
the magnetic fields. It means any realistic model for self- 
gravitating disk should consider the possible effects of the 
magnetic fields. Of course, our self-similar solutions are 
too simple to make any comparison with observations. 
But, we think, one may relax self-similarity assumption 
and solve the equations of the model numerically. In do- 
ing so, our self-similar solutions can greatly facilitate test- 
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ing and interpretation of results. Then, we can calculate 
the spectral energy distribution of such a self-gravitating 
magnetized disk. 
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Fig. 1. — Typical behaviours of the surface density (top, left), opening angle of the disk (top, right), ratios of the radial and the rotational 
velocities to the sound speed (middle), ratio of the thermal pressure to the magnetic pressure at the surface (bottom, left) and the Toomre 
parameters corresponding to a = 0.1, rjo = 0.1, = —0.8 and the other input parameters are indicated on the plots. Magnetic Toomre 
parameter Qm is denoted by the dotted lines. These plots show the effects of the mass accretion rate rh on the structure of the disk. 




Fig. 2. — The same as Figure 1, but for rh = 0.5, a = 0.1, and fd^ = —0.8. These plots show the effects of magnetic diffusivity coefficient 
770 on the solutions. Also, the magnetic Toomre parameter is denoted by dotted lines. 
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Fig. 3. — Top plots show the effects of the viscosity parameter a on the structure of the disk. Among the physical variables describing 
the disk, however, the ratio of the thermal pressure to the magnetic pressure at the surface of the disk /3 and the Toomre parameters are 
more sensitive to the variations of the viscosity coefficient. We have similar behaviour when the parameter f3ip varies, i.e. the ratio /3 and the 
Toomre parameters Q and Qm are more sensitive to the variations of keeping the rest of the input parameters fixed (bottom plots). 



